NASA TECHNICAL NOTE 


l n 


i 



on 



| NASA IN M257 0 J_ 


H 

m 

o 


Loan copy: f 

r AFWL TECHNIC/ 
KIRTLAND AF 


o | 

b-» ! 
LJj 
UJ = 
-0 ■ 
03 S 
-J i 


09 

5 TO 


■<\RY 

>1 


3 U 


MINIMUM-MASS DESIGN OF FILAMENTARY 
COMPOSITE PANELS UNDER COMBINED 
LOADS: DESIGN PROCEDURE BASED 
ON SIMPLIFIED BUCKLING EQUATIONS 


W. Jefferson Stroud and Nancy Agranoff 

Langley Research Center 
Hampton, Va. 23665 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 0. C. • OCTOBER 1976 



TECH LIBRARY KAFB, NM 


0133187 


1. Report No. 

NASA TN D-8257 

4. Title and Subtitle 


2. Government Accession No. 


3. Recipient's Catalog No. 


5. Report Date 


MINIMUM-MASS DESIGN OF FILAMENTARY COMPOSITE 
PANELS UNDER COMBINED LOADS: DESIGN PROCEDURE 
BASED ON SIMPLIFIED BUCKLING EQUATIONS 

7. Author(s( 

W. Jefferson Stroud and Nancy Agranoff 

9. Performing Organization Name and Address 


O ctob er 1978 

6. Performing Organization Code 


8. Performing Organization Report No. 

L- 10770 

10. Work Unit No. 

505-02-42-01 


NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Note 


National Aeronautics and Space Administration 
Washington, DC 20546 

15. Supplementary Notes 


16. Abstract 


14. Sponsoring Agency Code 


An analytical procedure is presented for designing hat-stiffened and corrugated panels 
made of composite material and subjected to longitudinal (in the direction of the stiffeners) 
compression and shear loadings. The procedure is based on nonlinear mathematical pro- 
graming techniques and a simplified set of buckling equations. Design requirements con- 
sidered are buckling, strength, and extensional and shear stiffness. Studies show the effects 
of specified thickness, variation of cross-section dimensions, stiffness requirements, local- 
buckling boundary conditions, and the effect of combined compression and shear loadings. 


17. Key Words (Suggested by Author(s)) 
Stiffened panels 
Panel design 
Structural synthesis 
Structural design 
Minimum-mass design 
Composite panels 

19. Security Classif. (of this report) 

Unclassified 



18. Distribution Statement 


Composite structure 
Structural efficiency 

Unclassified — Unlimited 

Subject Category 39 

1 20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 


51 

$4.25 


For sale by the National Technical Information Service, Springfield, Virginia 221 61 



CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 2 

LOADINGS AND CONFIGURATIONS CONSIDERED 5 

Loadings 5 

Configurations 6 

ANALYSIS -DESIGN PROCEDURE 7 

Stress Analysis 8 

Longitudinal compression 8 

Shear 8 

Stresses and strains in each lamina 9 

Strength Failure Criteria 9 

Maximum strain 9 

Stress interaction 10 

Buckling Analysis 10 

Overall buckling 10 

Local buckling 12 

Sizing Technique 13 

DESIGN STUDIES CARRIED OUT WITH ANALYSIS -DESIGN PROCEDURE 14 

Structural Efficiency Diagrams and Scaling 14 

Hat-Stiffened Panel, Longitudinal Compressive Loading N x 15 

Effect of discrete values of thickness tgt 45°) 15 

Ny 

Variation of dimensions as function of — for t2 = 0.56 mm (0.022 in.) . . 17 

Effect of changes in width of element 3 20 

Effect of changes in width of element 4 21 

Effect of local -buckling boundary conditions 23 

Stiffness considerations 25 

Hat-Stiffened Panel, Combined Shear and Compression 27 

Corrugated Panel, Longitudinal Compressive Loading 30 

Corrugated Panel, Combined Shear and Compression 31 

Hat- Stiffened and Corrugated Panels, Shear Loading 32 

CONCLUDING REMARKS 34 

APPENDIX: STIFFNESS EXPRESSIONS APPEARING IN ANALYSIS -DESIGN 

PROCEDURE 36 

iii 



Page 

Basic Laminate Constitutive Relations 36 

Extensional Stiffness 37 

Shear Stiffness 37 

Smeared Stiffnesses for Overall Buckling 38 

Bending stiffnesses El and Dj 39 

Bending stiffness Dg 40 

Stiffness Dg 40 

REFERENCES 42 

TABLES 44 


iv 


■ i ii i ii 


IHI llll I I 


I I I 


J 


MINIMUM-MASS DESIGN OF FILAMENTARY COMPOSITE PANELS 
UNDER COMBINED LOADS: DESIGN PROCEDURE BASED 

ON SIMPLIFIED BUCKLING EQUATIONS 

W. Jefferson Stroud and Nancy Agranoff 
Langley Research Center 

SUMMARY 

An analytical procedure is presented for designing hat-stiffened and corrugated 
panels made of composite material and subjected to longitudinal (in the direction of the 
stiffeners) compression and shear loadings. The procedure is based on nonlinear mathe- 
matical programing techniques and a simplified set of buckling equations. Design require- 
ments considered are buckling, strength, and extensional and shear stiffness. Studies 
show the effects of specified thickness, variation of cross-section dimensions, stiffness 
requirements, local-buckling boundary conditions, and the effect of combined compression 
and shear loadings. 


INTRODUCTION 

Because composite materials offer the potential for substantial mass and cost 
savings, Langley Research Center is engaged in a broad program to study structural 
applications of filamentary composite materials. One objective of that program is to 
establish a mass and strength data base for generic structural components made of com- 
posite material. Such a data base, containing both analytical and experimental results, is 
necessary to develop confidence in composite structural design procedures. The generic 
component considered in this report is the stiffened panel: specifically, hat-stiffened and 
corrugated panels. 

Some of the previous work dealing with the analysis and design of stiffened composite 
panels is discussed in references 1 to 6. The present report presents an analytical pro- 
cedure for designing low-mass, hat-stiffened and corrugated panels which are subjected to 
longitudinal compression and shear loadings. Design studies carried out with the pro- 
cedure are also presented. 

The procedure is based on nonlinear mathematical programing techniques in which 
buckling, strength, and stiffness requirements as well as upper and lower bounds on cross- 
section dimensions are taken into account with inequality constraints. The objective is to 
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obtain minimum-mass designs. Overall -buckling loads are calculated by assuming 
smeared stiffnesses. Local-buckling loads are calculated by assuming simple -support 
boundary conditions along the lines of attachment of adjacent elements. Certain buckling 
modes (stiffener rolling, for example) are neglected in the buckling analyses used herein. 
Lamina stresses and strains are used in the strength calculations. 

In composite-panel designs, the interaction of design variables may be quite com- 
plex. To demonstrate typical design study options available with the procedure presented 
in this report, studies were performed on hat-stiffened and corrugated graphite- epoxy 
panels loaded in compression and shear. These studies show the effects of discrete 
thickness, variations of cross-section dimensions, stiffness requirements, local- buckling 
boundary conditions, and the effect of combined compression and shear loadings. Detailed 
study results are presented so that a designer may obtain a better understanding of the 
mass penalties associated with practical design constraints such as stiffener cap widths 
and stiffener spacing. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The calculations were made 
in U.S. Customary Units. 


A 

A 


surface area of panel, pL 


area enclosed by hat (See eq. (A25).) 


i] 


laminate stiffnesses defined in equations (Al) and (A2) 




value of A ij for element k 


panel width (See fig. 5.) 


bk 

D 


element lengths defined in figure 2 


depth of element 2 defined in figure 23 


D ii 


laminate stiffnesses defined in equations (Al) and (A3) 


D ii, 


value of Dy for element k 


2 



smeared orthotropic stiffnesses defined in equations (A10) to (A13) 


d 1» d 2> d 3’ 

D x> D y> D xy 


d 

E 

El 


ET 

ET k 

e 1> e 2 


G 

GT 

G 


12 


N„ 


^x,crit 


N, 


N- 


xy 

•^xy,crit 


N- 


x y t 


N t 


distance between stiffeners; 2b^ 

Young's modulus 

beam bending stiffness of one period of panel cross section; defined in 
equation (A19) 

longitudinal extensional stiffness of panel; defined in equation (A5) 

longitudinal extensional stiffness for element k, defined in equation (A4) 

Young's modulus of composite material in fiber direction and transverse to 
fiber direction, respectively 

shear modulus 

shear stiffness of panel; defined in equation (A6) 

shear stiffness of composite material in coordinate system defined by fiber 
direction 

panel length 

applied longitudinal compressive loading per unit width of panel (See fig. 1.) 

value of compressive load that causes buckling 

longitudinal compressive loading in element k (See fig. 6.) 

applied shear loading per unit width of panel (See fig. 1.) 

value of shear load that causes buckling 

shear loading in element k (See fig. 6.) 

applied transverse (y-direction) loading per unit width of panel 



p 


period of stiffened panel; + 2b 4 


Q 11 ’ Q 12^1 
Q 22’ Q 66 J 

t k (e) 

w 

W/A 

L 

w 

x,y 


z k 


a 


xy 


e x’ e y» y xy 


e l ,e 2’ y 12 
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P 

M X jMy 


lamina stiffnesses 

thickness of lamina in element k at fiber orientation angle 0 (See fig. 2.) 
mass of one period of stiffened panel 
mass index 

lateral deflection of panel 

coordinate directions; axes defined in figure 1 

vertical distance from neutral axis of element 3 to neutral axis of panel cross 
section 

vertical distance from neutral axis of element 3 to neutral axis of element k 
angle element 2 makes with the horizontal (See fig. 23.) 
shear buckling parameter 

average overall shear strain of panel, defined by equation (A7) 
laminate strains in X-Y coordinate system 
laminate strains for element k 

lamina strains in coordinate system defined by fiber direction 

shear buckling parameter defined in equations (20) and (25), and fiber 
orientation angle shown in figure 1 

Poisson's ratio 

Poisson's ratios associated with bending; defined by equation (A10) 
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^12 ,(i 21 Poisson's ratios of composite material in coordinate system defined by fiber 
direction 

p density 

° r l> cr 2 >T ^2 lamina stresses in coordinate system defined by fiber direction 

r shear stress 

Superscript: 

a allowable 

Subscript: 

k integer denoting element number 

LOADINGS AND CONFIGURATIONS CONSIDERED 

Loadings 

Two loadings are considered: a uniform, longitudinal compressive loading N x in 
the direction of the stiffeners and a uniform shear loading Nxy. These loadings are 
shown in figure 1. 



Figure 1.- Loadings considered. 
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Configurations 

The two panel cross-sectional configurations considered in this report are shown 
in figure 2. The cross-section geometry is defined in terms of the four dimensions bj, 
b 2 , bg, and b^. The quantities bj, b 2 , and bg are the widths of elements 1, 2, and 3, 
respectively. The quantity b 4 is the half -width of element 4. Because of symmetry, 
b 3 

b 4 is replaced by -y in the corrugated panel. 

The panels are assumed to be made of a filamentary composite material with ply 
orientation angle 6 measured from a line drawn in the direction of the stiffeners (the 
x-direction) as shown in figure 1. In the examples presented in this report, only two 
fiber orientations are considered: 0=0° and 0 = ±45°. The 45° plys are assumed to 
be balanced and symmetric and are, therefore, denoted ±45°. 


p 



Figure 2.- Configurations and design variables. Element identification 
numbers are indicated in figure 2(a). 
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(b) Corrugated panel. 


and mmsm m e=±45° 

MMMMM333EM e= o° 

(c) Filament orientation pattern. 

Figure 2.- Concluded. 

All elements contain ±45° material. Elements 3 and 4 can also contain 0° material. 
The configurations are defined so that the 0° material and the ±45° material have specific 
locations. In both configurations t 3 (0°) refers to the thickness of one of the two equal 
0° layers in element 3; t 4 (0°) is the thickness of a 0° layer sandwiched between ±45° 

skins in element 4. 

The hat-stiffened configuration is a uniformly thick corrugated sheet of ±45° material 
attached to an equally thick flat sheet of ±45° material with 0° material at specified loca- 
tions in elements 3 and 4. The corrugated configuration is a uniformly thick corrugated 
sheet of ±45° material with 0° material at specified locations in element 3. In both 
configurations, the thickness of a ±45° sheet is denoted t 2 (±45°) or simply t 2 . The 
thicknesses of the ±45° sheets are taken to be uniform to provide a design that is relatively 
easy to fabricate. 


ANALYSIS-DESIGN PROCEDURE 

Nonlinear mathematical programing techniques are used to adjust the element widths 
and thicknesses (design variables) to provide minimum-mass designs that meet the design 
requirements. The design requirements (constraints) considered are strength, buckling, 
stiffness, and upper and lower bounds on the cross-section dimensions. The objective of 
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this section is to present and discuss briefly many of the expressions needed to evaluate 
the strength and buckling constraints. Equations for the extensional and shear stiffnesses 
are given in the appendix together with detailed expressions for other quantities that are 
used in this report. 


Stress Analysis 

The longitudinal and transverse strains and the longitudinal stress resultants are 
calculated for an applied longitudinal compressive load N x . The transverse loading N y 
is assumed to be zero for the overall panel and for each element that makes up the panel. 
The shear -stress resultants and shear strains are then calculated for an applied shear 
loading Njjy. Finally, by use of coordinate-transformation and stress-strain relations, 
the strains and stresses are calculated in the directions parallel and perpendicular to the 
fiber direction in each lamina of each element. All elements are assumed to be specially 
orthotropic (ref. 7) so that the inplane problem in each element is given by equation (A2) 
in the appendix. 

The stress analysis equations are summarized in table I and are discussed in the 

text. 

Longitudinal compression . - All elements of the panel cross section are assumed to 

have the same longitudinal strain e x given by equation (1). An expression for the 

longitudinal extensional stiffness ET of the panel cross section is given in the appendix. 

Since the transverse loading is assumed to be zero, the transverse strain e v in the kth 

y k 

element is given by equation (2). The longitudinal loading N x in the kth element is 

1C 

proportional to the longitudinal extensional stiffness ET^ of that element and is given 
by equation (3). 

Shear. - In the case of a hat-stiffened panel, the applied shear-stress resultant 
is assumed to be distributed over the panel cross section in the following manner. Ele- 
ment 4 carries the entire applied shear stress. In the region of the hat stiffener, the 
applied shear stress divides into two parts. Element 1 carries one part of the load, and 
elements 2 and 3 carry the other part. The relative proportions into which Nxy divides 
are obtained by equating the shear displacement of element 1 to the sum of the shear 
displacements of elements 2 and 3 as shown in figure 3. The resulting internal shear- 
stress resultants are given by equations (4) to (6), and the corresponding shear strains 
are given by equation (7). For a corrugated panel, the stiffness Agg^ is zero and 

N xy 2 = N xy g = N xy- 

During the synthesis, the internal stress resultants N x ^ and given by 

equations (3) to (6), are compared with the local -buckling stress resultants discussed in 
a subsequent section. 
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Figure 3. - Equating shear displacement of 
element 1 with sum of shear displacements 
of elements 2 and 3 to obtain internal shear - 
stress distribution. 

Stresses and strains in each lamina .- With 9 defined as the angle between the fiber 
direction and the x-direction (fig. 1), the strains’ in each lamina in the 1-direction, which 
is parallel to the fiber direction, and in the 2 -direction, which is perpendicular to the fiber 
direction, are given by equations (8) to (10). The corresponding stresses in each lamina 
are given by equations (11) to (13). 

The stresses and strains given by equations (8) to (13) are used in the strength 
criteria discussed in the next section. Since each element contains fibers oriented in both 
the -9 and +0 directions, stresses and strains at both orientations should be examined 
for strength considerations if shear loads exist. 

Strength Failure Criteria 

Several strength criteria have been proposed for designing filamentary composite 
laminates. Two criteria, referred to here as the maximum -strain criterion and the 
interactive-stress criterion, were used in the studies included in this report. (In no case 
were both criteria used simultaneously.) Both approaches appear, for example, in refer- 
ence 8. The chosen failure criterion is applied to each lamina in each element. 

The strength failure criteria equations are summarized in table II and are discussed 
in the text. 

Maximu m strain. - In the maximum -strain criterion, upper and lower limits are 
placed on the permissible values of e^, and u PP er limit is used with 
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positive strains, and the lower limit with negative strains. In terms of inequality expres- 
sions, the maximum -strain criterion is given by equations (14) to (16) in which the super- 
script a means allowable, and signs and values of allowable strains must be selected to 
limit both positive and negative strains. 

Stress inte raction . - While each separate strain component is limited in the 
maximum -strain criterion, two of the stress components are linked in the stress- 
interaction criterion given here. Fiber failure is assumed to be prevented by equa- 
tion (17), and matrix failure is assumed to be prevented by equation (18). Again, the 
superscript a means allowable, and signs and values of allowable stresses must be 
selected to limit both positive and negative stresses. 

Buckling Analysis 

The approach taken here is a traditional one in which buckling is assumed to be 
prevented if several independent forms of buckling, including overall- and local-buckling 
modes, are postulated and prevented. Overall -buckling equations are based on smeared 
stiffnesses, and local -buckling equations are based on simple -support boundary conditions 
along the lines of attachment of adjacent elements. (This same approach was used, for 
example, in refs. 9 and 10 and more recently in ref. 3.) This approach ignores compati- 
bility of buckling mode shape, a factor that tends to make the approach conservative. On 
the other hand, the simplified analyses do not consider all the potential buckling modes, a 
factor that could make the approach unconservative. Unpublished studies with a more 
accurate buckling analysis indicate that in most cases the approach is sufficiently accurate 
to predict structural efficiencies and design trends. In some cases, such as those in 
which roll modes are active, the buckling analysis approach used here is not adequate. 

Both the overall- and local-buckling analyses assume that the plate structure is 
orthotropic with principal axes alined with the edges of the plate. Buckling caused by 
longitudinal compression and buckling caused by shear are treated separately. Combined 
loads are accounted for with a buckling-interaction formula. 

The buckling equations, which are based on work presented in references 8 and 11 
to 15, are summarized in table III and are discussed in the text. Certain detailed expres- 
sions, including equations for the smeared orthotropic stiffnesses, are presented in the 
appendix. 

Overall buckling. - The order of presentation is longitudinal compression, shear, and 
combined longitudinal compression and shear. 

Longitudinal compression: The stiffened panel is treated as a wide column with the 
loaded ends simply supported arid with the lateral edges free as shown in figure 4. A 
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buckling equation that takes into account transverse shear effects is presented as equa- 
tion (19) in table III. The transverse shearing-stiffness term is taken as 2 * b 2 * Agg^ 



Figure 4.- Wide-column or Euler buckling is assumed 
for overall buckling in axial compression. 


Shear: The panel is assumed to be long in the y -direction and simply supported 

along the edges x = 0 and x = L as shown in figure 5. The buckling equations are 
given as equations (20) to (22) in table III. 

Equations (21) and (22) are based on the assumption that the smeared bending stiff- 
ness in the x-direction is much larger than the smeared bending stiffness in the 

y-direction Dg, and on the assumption that the dimension a (fig. 5(a)) is larger than 
the dimension L. These assumptions permit a parameter /3 defined by 
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to be taken as zero in the derivation of equations (21) and (22). The equations correspond, 
therefore, to the infinitely long panel shown in figure 5(b). 



SS= Simple support 

(a) Finite shear model. (b) Infinitely long shear model. 

Figure 5. - Finite and infinitely long shear models for overall shear buckling. 


Combined longitudinal compression and shear: When both longitudinal compression 
and shear loadings are present, buckling is assumed to occur when the buckling-interaction 
formula given by equation (23) is satisfied. Based on studies of isotropic plates in refer- 
ence 14 and orthotropic plates in reference 15, this buckling-interaction formula appears 
to be adequate for the panels and loadings considered herein. 

Local buckling . - The second level subscript (indicating the kth element) is dropped 
in table III and elsewhere in this report except where it is necessary for clarity. 

Longitudinal compression: The structural model is shown in figure 6. The length 
L of each element is assumed to be much greater than its width b. With these simpli- 
fying assumptions the buckling equation, which is applied to each of the four elements of 
the cross section, is given by equation (24). 

Shear: Local shear buckling is treated in the same way that overall shear buckling 

was treated in an earlier section. The only difference is the stiffnesses that are used. 
Again, the structural model is shown in figure 6. The equations for local shear buckling 
are given by equations (25) to (27) in table III. 

Equations (26) and (27) are based on the assumption that L » b (fig. 6). This 
assumption permits a parameter / 8 defined by 



1/4 


P 


b/ P ll 

L V D 22 


to be taken as zero in the derivation of equations (26) and (27). 



Combined longitudinal compression and shear: The buckling-interaction formula 

given by equation (23) is used for local as well as overall buckling. The adequacy of that 
formula is examined and verified for orthotropic plates in reference 15. 

Sizing Technique 

The design problem is formulated as a nonlinear mathematical programing problem. 
The design variables are the element widths and thicknesses that define the cross section. 
The constraints are the strength, buckling, and stiffness design requirements. Since the 
panel length L is fixed, the objective function (the quantity that is minimized) is the 
panel mass per unit width. Equations for the strength and buckling constraints have been 
presented in the preceding sections. Equations for the stiffnesses are given in the 
appendix. 
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The specific design constraints considered herein are overall buckling, local 
buckling of each element, strength of each lamina in each element, extensional stiffness, 
and shear stiffness. Upper and lower bounds are also placed on various dimensions. 

The computer program used to perform the optimization in the studies presented 
herein is a general-purpose parameter optimizer denoted AESOP (refs. 16 to 18) which 
accounts for constraints using an exterior penalty function approach. The AESOP code 
contains several optimization algorithms that can be selected in various combinations by 
the user. The approach used here was a combination called "creeper" and "pattern." 
These algorithms are described in references 16 to 18. 

DESIGN STUDIES CARRIED OUT WITH ANALYSIS-DESIGN PROCEDURE 

Studies were carried out to examine detailed design trends, sensitivity of structural 
efficiency to changes in certain dimensions, effect of extensional and shear -stiffness 
design requirements, effect of local buckling boundary conditions, and the effect of com- 
bined compression and shear loadings. The material properties that are used for all 
studies are given in tables IV and V. The upper and lower limits placed on the cross- 
section dimensions are given in table VI. 

Structural Efficiency Diagrams and Scaling 

The structural design data developed in this report are presented in the form of 
structural efficiency diagrams in which the mass index of minimum -mass panels 

/n x NjJ) l 

is given as a function of the loading index or -jf-J, also called the structural index. 
Presenting design data in this manner minimizes size effects and broadens the applica- 
bility of the data. If there were no upper or lower limits on the cross-sectional dimen- 
sions and if the thickness could vary continuously, then size effects would be eliminated. 
Size effects are eliminated because of dimensional similarity or scaling principles that 
apply to composite panels in the same way that they have been applied to conventional 
metal panels. (See, for example, ref. 19.) 

The scaling principle for panels with loadings considered herein can be summarized 
as follows. If two panels have the same proportions - that is, if all the dimensions 
(including thicknesses and length) of one panel are multiplied by a single scale factor to 
obtain the dimensions of the other panel - and if both panels are subjected to loads which 
give the same value of the loading index, then 

(1) the stresses in the corresponding laminas are the same for both panels; 

(2) the ratios of the element stress resultants to the element buckling-stress 

resultants are the same for both panels; 
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(3) the ratio of the overall stress resultant to the overall buckling-stress resultant 

is the same for both panels; and 

(4) the mass index is the same for both panels. 

These four statements confirm that for a given structural concept (including the material 
system), minimum-mass designs fall on a single structural efficiency curve and that 
dimensions of minimum-mass panels scale linearly with the load (not the load index). 

The results are presented in terms of dimensional quantities. In some cases, for 
example, t2 is used as a parameter. In all cases the data are for panels with 
L = 0.76 m (30 in.). It is understood that scaling techniques can be applied to the results 

to make them applicable to other panel sizes and loads. For example, curves for which 

to 

t 2 is constant can be interpreted as curves for which is constant. 

Hat-Stiffened Panel, Longitudinal Compressive Loading N x 
Effect of discrete values of thickness t 2 (±45°).- Minimum values of the mass index 
are given in figure 7 as a function of compressive loading index for various 

T 1 ^ i_i 

values of t2(±45°). It is assumed that the thickness of 0 material can vary continuously. 
The values chosen for t 2 represent four-ply and eight-ply laminates made up of plies 
of either 0.051 mm (0.0020 in.) or 0.140 mm (0.0055 in.) thickness. Laminates that are 
multiples of four plys can be made balanced and symmetric; balanced and symmetric 
laminates are required in all the studies presented in this report. The curves give an 
indication of the range over which a particular value of t 2 is lightest and the penalty 
involved in selecting another thickness. 

The dashed curve indicates designs for which t 2 (±45°) is allowed to vary contin- 
uously. That curve is referred to as the t 2 free curve. It is tangent to the various 
curves for which t 2 is constant. The designs associated with the tgfree curve are 
referred to as the t 2 free designs. 

The requirement that elements 1 and 3 have a width not less than 2.03 cm (0.80 in.) 

begins to have an effect for loadings less than = 276 kPa (40 lbf/in 2 ) where the 
tgfree curve begins to deviate from a straight line. 

All the examples presented herein use the interactive -stress criterion given by 

In- 
equations (17) and (18). For loadings greater than — - = 3.45 MPa (500 lbf/in 2 ), strength 

ij 

considerations begin to affect the t 2 free designs and the designs for which t 2 = 1.12 mm 
(0.044 in.). For this hat-stiffened configuration and for the material properties (table IV) 
used for these studies, the maximum strain requirement (eqs. (14) to (16)) produces 
slightly heavier designs than the interactive-stress requirement. For example, at 
N 

— = 6.89 MPa (1000 lbf/in 2 ) where the difference is largest, the maximum strain designs 
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are about 4 percent heavier than the interactive-stress designs. This is the only study 
reported in this paper in which the maximum -strain criterion is used. Strength require- 
ments do not affect the designs with the other values of t 2 shown in figure 7. 



Figure 7.- Structural efficiency of graphite -epoxy, hat-stiffened compres- 
sion panels for various prescribed values of thickness t£ and for t 2 
varying continuously; L = 0.76 m (30 in.); thickness t 2 is shown in 
mm (in.). 

Consider the loading ^ = 689 kPa (100 lbf/in 2 ). The lightest design has t 2 
equal to about 0.41 mm (0.016 in.). Manufacturing requirements may, however, make it 
more desirable to fabricate a design with t 2 equal to 0.56 mm (0.022 in.). The cross 
sections for the minimum-mass designs having t 2 equal to 0.41 mm and 0.56 mm are 
presented in figure 8. The design for which t 2 = 0.56 mm is larger and deeper than 
the design for which t 2 = 0.41 mm. (Both designs are penalized slightly by the 
bg 5 2.03 cm (0.80 in.) requirement.) 
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L 

Inches 0 

L 

Centimeters 0 



tg = 0.41 mm (0.016 in.) 

tg = 0.56 mm (0.022 in.) 


Figure 8.- Cross sections for two hat-stiffened panels 

designed for ^ = 689 kPa (100 lbf/in2); L = 0.76 m 
(30 in.). 

The significant point here is that even though t 2 is increased from 0.41 mm 
(0.016 in.) to 0.56 mm (0.022 in.) - a 37-percent increase - the design for which 
t£ = 0.56 mm is only 4-percent heavier than the design for which t 2 = 0.41 mm. In this 
case, as well as in many other cases, the manufacturing requirement on the thickness t 2 
can be accommodated with very little penalty in structural efficiency. It should be pointed 
out, however, that although there is very little difference in the mass of the two designs, 
the extensional and the shearing stiffnesses of the one design differ substantially from the 
corresponding stiffnesses in the other design. This fact is discussed in more detail in a 
subsequent section. 

N v 

Variation of dimensions as function of for t 2 = 0.56 mm (0.022 in.). - Cross- 

section dimensions of the designs in figure 7 for which t 2 = 0.56 mm (0.022 in.) are 

N 

presented as a function of the loading index -=^ i n figure 9. These dimensions are for 
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Figure 9.- Variation of cross-sectional dimensions as function of 
loading for minimum-mass, hat -stiffened, graphite -epoxy 
compression panels; t 2 = 0.56 mm (0.022 in.); L = 0.76 m 


a panel with a length L of 0.762 m (30 in.) and for the material properties of table IV. 
The loading ranges over which various buckling constraints are active are given at the top 
of figure 9. Element 3 is not buckling critical over the loading range shown, in figure 9. 

The width b 3 of element 3 is shown in figure 9 as being equal to its minimum 

N 

value (2.03 cm (0.80 in.)) until — - is equal to about 4830 kPa (700 lbf/in 2 ). Actually, 

L N 

b 3 has a weak tendency to increase to a larger value over the range = 2410 to 
3450 kPa (350 to 500 lbf/in 2 ). This increase in bg is accompanied by a decrease in 
tg(0°) so that at a given loading the total longitudinal extensional stiffness (ETg • bg) of 
element 3 remains almost constant regardless of the relative values of bg and tg(0°). 
The simple variation of bg and tg(0°) is shown in figure 9 for clarity. The weak 
tendency of bg to exceed its minimum value in the above loading range is discussed in 
the next section. 

Cross-section shapes of the panels having tg = 0.56 mm (0.022 in.) and designed to 
N v 

support loadings of — = 68.9, 689, and 6890 kPa (10, 100, and 1000 lbf/in 2 ) are pre- 

L N 

sented in figure 10. The designs for — — = 68.9 and 6890 kPa are shown merely to 

l_i 

indicate the design trends for a fixed value of tg. Other values of tg provide substan- 
tially lighter designs at these two loadings. 

Inches 0 1 i j 

0 2 4 6 8 

Centimeters I ! I ! i i I i j 




68.9 kPa (10 lbf/in 2 ). 



(c) = 6890 kPa (1000 lbf (in 2 ). 

Figure 10.- Cross-sectional shapes of hat- stiffened, graphite -epoxy compression 
panels having t 2 = 0.56 mm (0.022 in.) and L = 0.76 m (30 in.) for three 
values of axial compressive load. 
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Effect of change s in width of element 3. - In the previous section it was pointed out 
that although the width b„ of element 3 is shown in figure 9 as being equal to the mini- 

mum value (2.03 cm (0.80 in.)) until ~ is equal to about 4830 kPa (700 lbf/in I 2 ), there 
is actually a weak tendency for bg to be larger than the minimum value in the region 

= 2410 to 3450 kPa (350 to 500 lbf/in 2 ). That effect is shown in figure 11 for a loading 

k n 

of = 2760 kPa (400 lbf/in 2 ). In figure 11 stiffened panels are completely redesigned 
Li 

for various prescribed values of bg; that is, each point shown by a circular symbol 
represents a minimum-mass design for that value of bg. For these calculations, 
t£ = 0.56 mm (0.022 in.). The minimum value of the mass index occurs at bg 

equal to about 2.54 cm (1.0 in.). For this loading, the penalty for choosing bg = 2.03 cm 

N x 

is about 0.25 percent. It is apparent that for this value of — , the width bg can vary 
over a large range with very little change in structural efficiency. The width bg is a 
weak design variable for other loadings also. This fact can be exploited to accommodate 
manufacturing or other design considerations. 



i i i i i .. . i 

2.0 2.5 3.0 3.5 4.0 4.5 cm 

I I L I i __ l 

.8 1.0 1.2 1.4 1.6 1.8 in. 

b 3 

Figure 11.- Percent change in structural efficiency as function of 
dimension bg for minimum-mass, graphite -epoxy, hat-stiffened 
compression panels; tg = 0.56 mm (0.022 in.); L = 0.76 m 

(30 in.); = 2760 kPa (400 Ibf/in 2 ). 
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Effect of changes i n width of element 4. - The structural efficiency is more sensitive 
to changes in than to changes in bg. An indication of that sensitivity is shown in 
figure 12 in which the percent change in mass index is presented as a function of width d 

N-* 

(twice b^ of the flat portion between stiffeners for two values of -jp: 689 and 4830 kPa 
(100 and 700 lbf/in 2 ). The circular and square symbols indicate minimum-mass designs 
for the 689 and 4830 kPa cases, respectively. For these results, the thickness t£ is 
allowed to vary freely. 
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Percent change in 
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Figure 12.- Effect on structural efficiency of increasing stiffener 
spacing d of hat -stiffened, graphite-epoxy compression panels; 
t2 varies continuously; L = 0.76 m (30 in.). 

N 

First, consider the case for which = 689 kPa (100 lbf/in 2 ). With no constraint 
placed on d, the lightest design occurs for d equal to about 4.8 cm (1.9 in.). As d 
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is increased, the mass index — j — increases. When d is equal to about 7.9 cm 

X j n 

(3.1 in.), the lightest design changes from one with 0 U material in elements 3 and 4 to one 
that has 0° material in element 3 only ^(0°) = 0^. The two concepts, one with and the 
other without 0° material in element 4, do not evolve from one another. They are com- 
pletely independent. The concept with no 0° material in element 4 is similar to those for 
low loadings shown in figure 9 in that bj is equal to b^ As d is increased from 
about 4.8 cm to 12.7 cm (5.0 in.), the mass index increases by about 42 percent. 

For the case in which = 4830 kPa (700 lbf/in 2 ), the efficiency of the concept for 
which t 4 ( 0 °) =0 is restricted because of the upper bound placed on t 3 (0°) and because 
of strength considerations. For these reasons, the concept with 0° material in both 
elements 3 and 4 is lighter than the concept with 0° material in only element 3 for the 
entire range of d considered. The lightest design occurs for d equal to about 8. 1 cm 
(3.2 in.). As d is increased to 12.7 cm (5.0 in.), the mass index increases by about 
18 percent. 

Another way of studying the effect of increasing the distance d between stiffeners 
is with a structural efficiency diagram such as that shown in figure 13. The objective 
here is to show the effects of a constraint defined by d £ 12.7 cm (5.0 in.). Since there 
is no tendency for d to be greater than 12.7 cm for the cases considered, the constraint 
becomes simply d = 12.7 cm. 

The bottom solid curve in figure 13 gives the structural efficiency of hat-stiffened 
designs with no constraint on d. It is the t 2 free curve of figure 8. The two dashed 
curves that are tangent to the t 2 free curve are portions of the corresponding curves 
of figure 8. Like the solid curve, these two dashed curves have no constraint on d and 
are used for reference. For the other curves, d is required to be 12.7 cm (5.0 in.). 

The other solid curve in figure 13 forms the lower bound for the designs having 
d = 12.7 cm (5.0 in.). For this curve, t2 is allowed to vary freely. For loadings less 

than about = 3450 kPa (500 lbf/in 2 ), the designs for which t^O 0 ) = 0 are lighter than 
the designs for which t^O 0 ) + 0. However, because of the upper bound on t3(0°) and 

because of strength considerations, the designs for which t^O 0 ) * 0 are lighter for 

N 

loadings greater than about = 3450 kPa (500 lbf/in 2 ). 

The penalty in structural efficiency caused by increasing the distance d between 

N x 

stiffeners to 12.7 cm (5.0 in.) depends upon the loading level. At — = 68.9 kPa 

L 

9 N X 

(10 lbf/in* 1 ), the increase in the mass index is about 77 percent; at — = 6890 kPa 
(1000 lbf/in 2 ), the increase is about 14 percent. 
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-p- . kPa 

Figure 13.- Effect of d = 12.7 cm (5 in.) stiffener spacing requirement on 
structural efficiency of hat -stiffened, graphite -epoxy compression panels. 
Thicknesses are shown in mm (in.); L = 0.76 m (30 in.). 

Effect of local -buckling boundary conditions. - As explained earlier, the analysis for 
local buckling is based on the assumption that each element is simply supported along its 
line of attachment to the adjacent elements. This assumption is believed to be conserva- 
tive in most cases. In particular, since element 3 never experiences local buckling in the 
designs generated in this report, it is reasonable to assume that element 3 provides rota- 
tional support to element 2 and that a more sophisticated analysis-design procedure could 
exploit that extra support. 

To assess the effect of local -buckling boundary conditions and to provide information 
on the potential benefits available with a more sophisticated analysis, additional designs 
were generated under the assumption that the local-buckling boundary conditions for 
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elements 1 and 2 were clamped rather than simply supported. The results are presented 
in figure 14. The solid curves represent designs obtained with the clamped boundary 
condition assumption, and the dashed curves are the original designs with the simple- 
support boundary conditions. 
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Figure 14. - Effect of local -buckling boundary conditions on structural efficiency 
of hat- stiffened, graphite -epoxy compression panels; L = 0.76 m (30 in.). 
Thicknesses are shown in mm (in.). 


The benefits derivable from the clamped boundary condition depend upon the curves 
being compared and the loadings at which the curves are compared. In the loading region 
where the designs with t£ = 0.56 mm (0.022 in.) are the most efficient designs, clamping 
provides a 10-percent reduction in the mass index. These results suggest that a more 
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sophisticated buckling analysis, such as those described in references 20 and 21, would 
provide designs that are about 10 percent more efficient than those presented in this 
report. 

Stiffness considerations. - It was shown in reference 5 that the longitudinal exten- 
sional stiffness of minimum-mass graphite -epoxy compression panels can be less than 
that of minimum-mass aluminum panels designed to carry the same load. This fact is 
demonstrated in figure 15 in which the longitudinal extensional stiffness ET is presented 

N x 

as a function of the loading index -=— for a portion of the designs shown in figure 7. 


N x jbf 



Figure 15.- Longitudinal extensional stiffness as function of loading for 

minimum-mass, graphite -epoxy and aluminum hat-stiffened compression 
panels; L = 0. 76 m (30 in.). Thicknesses are shown in mm (in.). 


The solid lines represent graphite -epoxy panels with four values of t 2 , and the dashed 
line represents aluminum panels. The discontinuities in the extensional stiffness of the 
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graphite -epoxy designs occur where the most efficient designs shift from one value of t 2 
to another. This relationship can be seen by comparing figure 15 with figure 7. The main 
reason for the "sawtooth” effect in figure 15 is the variation of 0° material in the panels. 

If stiffness is an important design consideration, that fact can be taken into account by 
including stiffness constraints in the synthesis procedure. 

Aluminum wing panels on commercial jet transports may have a much higher longi- 
tudinal extensional stiffness than that shown in figure 15. (See ref. 6.) Commercial air- 
craft may also have large shear-stiffness requirements as well as stiffener spacing 
requirements. To study the effect of these additional design requirements, minimum- 
mass panels were designed for various combinations of extensional and shear stiffness. 
(Expressions for both stiffnesses are given in the appendix.) All panels were designed to 

carry a load of ^ = 689 kPa (100 lbf/in^), and all panels have a distance d between 

J-J 

stiffeners of 12.7 cm (5.0 in.) and a length L of 0.76 m (30 in.). The thickness t£ is 
allowed to vary freely. 

W /A 

The results are presented in figure 16 in which the mass index — ^ — of the 

L 

minimum-mass panels is shown as a function of the design requirements on shear stiff- 
ness and longitudinal extensional stiffness. The dashed lines are lines of constant value 
of the mass index. The circular symbol represents the design which has no stiffness 
requirement. The solid line above the symbol represents the locus of minimum-mass 
designs that meet an extensional- stiffness. requirement, and the solid line below and to the 
right of the solid symbol represents the locus of minimum-mass designs that meet a 
shear -stiffness requirement. Points within the region containing the dashed lines repre- 
sent designs with combined extensional- and shear -stiffness requirements. The discon- 
tinuity in the line defining the extensional-stiffness designs is caused by a change in 
design concept similar to that shown in figure 12. 

Two observations follow from the results presented in figure 16. First, starting 
at the symbol, increasing either stiffness with a minimum increase in mass (that is, 
moving along the solid line) causes the other stiffness to decrease. Except for the lower 
boundary, this characteristic exists for all designs. This characteristic exists primarily 
because the extensional stiffness is increased by adding 0° material and the shear stiff- 
ness is increased with ±45° material. Second, the penalty in structural efficiency 
associated with increasing the shear stiffness is much greater than the penalty associated 
with increasing the extensional stiffness. 
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Figure 16.- Panel structural efficiency as function of longitudinal-extensional- 
stiffness and shear-stiffness requirements for minimum-mass, hat -stiffened, 
graphite -epoxy compression panels designed to support loading of 

-jj = 689 kPa (100 lbf/in 3 ). Thickness t 2 varies continuously. Mass 
index — given in kg/m 3 (lbm/in 3 ). 

Hat-Stiffened Panel, Combined Shear and Compression 
The effect of a shear loading on the structural efficiency of hat-stiffened compres- 
sion panels is shown in figures 17 and 18. The mass index i s presented as a 

„ . N x Nxv ^ Nyv 

function of for several values of in figure 17 and for several values of ■ ■■■ 

. L N x 

in figure 18. For each curve, t 2 is allowed to vary continuously, so. the cusps associated 

with constant values of t 2 do not appear. As is to be expected, the larger the ratio of 
Nxy N x 

— jj- to -j-, the larger the effect on the mass index. 
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Figure 17.- Effect of shear loading on structural efficiency of hat -stiffened, 
graphite -epoxy compression panels; L = 0.76 m (30 in.); t2 varies 
continuously. 
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Corrugated Panel, Longitudinal Compressive Loading 

To study the effect of discrete thickness, corrugated compression panels were 
designed with various thicknesses of ±45° material. The results, which are presented in 
figure 19, show the same general trends as the hat-stiffened panels (fig. 7). The two 
main differences are: first, the values of the mass index are lower for the corrugated 
panel than the hat-stiffened panel; second, the most efficient loading range for a given 
value of t 2 is lower for a corrugated panel than for a hat- stiffened panel. 





-pkPa 

Figure 19.- Structural efficiency of corrugated, graphite -epoxy compression 
panels for various values of thickness t 2 ; L = 0.76 m (30 in.). 
Thicknesses are shown in mm (in.). 
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Corrugated Panel, Combined Shear and Compression 

The effect of a shear loading on the structural efficiency of corrugated compression 
panels is shown in figures 20 and 21. The curves are similar to those presented in 
figures 17 and 18. The main difference is that the mass index is lower for the corrugated 
panels than for the hat- stiffened panels. 


N y Jbf. 

— > 3 



-p-. kPa 

Figure 20.- Effect of shear loading on structural efficiency of corrugated, 
graphite -epoxy compression panels; L = 0.76 m (30 in.); t£ varies 
continuously. 
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Figure 21.- Structural efficiency of graphite -epoxy corrugated panels 
subjected to combinations of longitudinal compression and shear 
loadings; L = 0.76 m (30 in.); t£ varies continuously. 


Hat-Stiffened and Corrugated Panels, Shear Loading 

The structural efficiency of hat-stiffened and corrugated panels loaded only by a 
shear loading is shown in figure 22. The solid curves are for graphite -epoxy panels. 

Over the range of loading considered, the graphite -epoxy corrugated panels are 5 to 
20 percent lighter than the graphite -epoxy hat-stiffened panels. 

The dashed curve is for aluminum corrugated panels having uniform thickness and 
designed using the same procedure as that used for the composite panels. The properties 

assumed for the aluminum material are given in table V. In the lower loading range 
N 

— — < 345 kPa (50 lbf/in^), where the panels are buckling critical but not strength critical, 
L 

the aluminum corrugated panels are about 45 percent heavier than the graphite -epoxy 
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corrugated panels. For loadings greater than — - — = 345 kPa, material strength con- 

Li 

siderations become dominant for the aluminum panels, and these panels become consider- 
ably heavier than the graphite -epoxy panels. 



N xv 

, kPa 

Figure 22.- Effect of shear loading on structural efficiency of graphite -epoxy 
and aluminum corrugated panels and graphite -epoxy hat-stiffened panels; 
L = 0.76 m (30 in.); t 2 varies continuously. 


The symbol in figure 22 represents the theoretical structural efficiency of an 

aluminum corrugated shear panel given in reference 22. The structural efficiency of that 

panel is shown by the intersection of the dashed curves (overall- and local -buckling 

curves) at -^ = 0.05 ksi in figure 9 of reference 22. (The quantity is the same as 

N xv\ h2 \ h 

— The structural efficiency is influenced by a 12 -percent increase in mass caused by 
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doublers at the edges x = 0 and x = L and by edge constraints that provide more 
support than simple -support boundary conditions. That panel is about 3 percent heavier 
than the aluminum corrugated panel designed with the procedure presented in this report. 

CONCLUDING REMARKS 

An engineering design procedure for sizing composite hat-stiffened and corrugated 
panels subjected to longitudinal compression and shear loadings has been developed and 
exercised. The procedure is based on nonlinear mathematical programing techniques. 
Expressions are given for the design constraints used in the procedure. These constraints 
are: overall-buckling loads, local-buckling loads, lamina strength, longitudinal exten- 
sional stiffness, and shear stiffness. 

The buckling analysis is simplified in that only certain specific buckling modes are 
considered. Because explicit eigenvalue solutions are available for these buckling modes, 
eigenvalue analyses are not required. Design studies covering a broad range of design 
parameters can, therefore, be carried out with a relatively small amount of computer 
time. Unpublished studies with a more accurate buckling analysis indicate that in most 
cases the approach is sufficiently accurate to predict structural efficiencies and design 
trends. In some cases, such as those in which roll modes are active, the buckling analysis 
approach used here is not adequate. 

In order to obtain a clearer understanding of both the analysis-design procedure and 
the composite panel design results that could be obtained by use of that procedure, several 
design studies were carried out. Buckling and strength design requirements are included 
in all the studies. In one study the effect of stiffness requirements is explored. Struc- 
tural efficiency diagrams are used to present the mass for panels designed to support 
longitudinal compressive loads, shear loads, and the increase in mass required for 
combined shear and compressive loadings. The studies indicate that panel mass is not 
extremely sensitive to panel dimensions; that is, panels with dimensions moderately 
different from those of the lightest panels may be only slightly heavier than the lightest 
panels. This fact can be exploited to accommodate manufacturing and other design con- 
siderations. It was also determined that scaling techniques in common use for sizing 
conventional metal panels apply equally well for composite panels where the number of 
design variables is large. 

Other observations, based on limited studies, are: 

1. There is a greater mass penalty associated with increasing the shear stiffness of 
hat-stiffened panels than with increasing the extensional stiffness. 
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2. The mass of hat- stiffened panels is affected more by changes in the distance 
between stiffeners than in the width of the hat cap. 

3. As has been found in studies for compression panels, composite shear panels 
offer substantial mass savings over metal shear panels. 

The design procedure presented in this report should be used in conjunction with 
engineering- judgment and other analyses to insure the adequacy of the final design. Other 
design requirements arising from manufacturing considerations, for example, could also 
have an impact on the design. To the fullest extent possible, all quantifiable requirements 
should be considered early in the synthesis. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 23, 1976 
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APPENDIX 


STIFFNESS EXPRESSIONS APPEARING IN 
ANALYSIS-DESIGN PROCEDURE 

Basic Laminate Constitutive Relations 
The plate constitutive equation can be written as 
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where Aij, Bjj , and Djj are laminate stiffnesses that are calculated from basic mate- 
rial properties by using standard expressions found, for example, in reference 8. The 
quantities M x , My, and M X y are plate moments, and k x , ky, and kxy are the asso- 
ciated curvatures. In all cases considered in this report, the ply-stacking sequence is 
balanced and symmetric so that the Ajg, A 2 g, and Bjj terms are zero. Also, the 
D-jg and D 2 g terms are assumed to be zero. For these types of laminates, which are 
denoted "specially orthotropic," the set of equations represented by equation (Al) uncouple 
and reduce to 
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In subsequent equations, the notations and Dij^ denote the values of Aij and 

Djj, respectively, for element k. 

Extensional Stiffness 

If the transverse load Ny is zero, the longitudinal extensional stiffness ET k of 
element k is given by 

ET k = A n - ^ 12k ^ (A4) 

R k A 22k 

If it is assumed that all elements of the panel cross section have the same longitudinal 
strain e x , the longitudinal extensional stiffness ET is given by 


ETjbi + 2ET 2 b2 + ET 3 b 3 + 2ET 4 b 4 
ET = bj + 2b 4 

For the corrugated panel, ETj is zero. 

Shear Stiffness 

The shear stiffness GT of the panel is defined as 


(A5) 


N 

GT = - 


'xy 


(A6) 


xy 


in which is the average shear angle given by 

_ 2 rxys b 2 + rxy 3 b 3 + 2y xy 4 b 4 ( A7 ) 

7 xy bj + 2b 4 

With equations (4) to (7) used to calculate the individual shear angles, the shear stiffness 
GT is 


GT = 


bj + 2b 4 


2b 2 b 3 


2b 4 
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(A8) 
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Equation (A8) can be used for either hat-stiffened or corrugated panels. For hat-stiffened 
panels only, the following equation is equivalent to equation (A8) and provides more insight 
into the effect of the hat on the stiffness GT: 
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Smeared Stiffnesses for Overall Buckling 

The overall -buckling equations given in table III are based on the assumption that 
the stiffened panel is an orthotropic plate with smeared stiffness properties. The smeared 
stiffnesses are defined in the following differential equation for lateral deflections of an 
orthotropic plate loaded as shown in figure 1: 
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( A 10) 


The stiffnesses that appear in equation (A10) are related to the stiffnesses that appear in 
equations (20) to (23) by 
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The objective of this section is to provide expressions for El, D-^, D 2 , and Dg. 
The approach is the same as that used in reference 23. 

Calculations presented in reference 23 show that for metal corrugation- stiffened 
panels, the product can be neglected when compared with 1.0. It is assumed that 

the product can also be neglected for composite panels of the type considered in this 
report. 
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Bending stiffnesses El and Dj_.- The panel bending stiffness D x is simply the 

El per unit width of a beam having the cross-sectional shape of the panel. Since p x /Uy 
is small compared with 1.0, Dj is equal to D x . It is first necessary to locate the neu- 
tral axis of the cross section. The distance from the center of element 3 to the neutral 
axis of element k is denoted z^. These distances are shown in figure 23 and are sum- 
marized below. 



Figure 23.- Additional dimensions for hat- stiffened panel. 


Hat-Stiffened Panel 

Corrugated Panel 


z 1 = D + t 4 (0°) + t 2 (±45°) 

z 1 = D 

(A14) 

z 2 = D/2 

Z 2 = D/2 

(A15) 
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(A16) 

Z 4 = D + l[t 4 (0°) + t 2 (±45°)] 

Z 4 = D 

(A17) 


The distance z from the center of element 3 to the neutral axis of the panel cross sec- 
tion is then given by 

_ = Z ET k b k fk (A18) 

Z ET k b k 
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in which ET k is given by equation (A4). In equation (A18), proper account must be taken 
of the fact that elements 2 and 4 appear twice in one period, and the summations are over 
one period. The bending stiffness El of one period of the panel cross section is 

(A19) 

The bending stiffness of each element about its own centroid is neglected in equation (A19). 
The error introduced by this assumption is only about 1 percent. 

The smeared bending stiffnesses D x and D-^ are the El per unit width of the 

panel 


D X = D x = f (A20) 

in which p is the period b^ + 2b^. 

Bending stiffness D2.- The approach presented in reference 23 is again followed. 
The smeared bending stiffness in the y-direction is given by 
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(A21) 


(A22) 


(A23) 


Angle a is the angle that element 2 makes with the horizontal as shown in figure 23. 

Stiffness D3.- For the hat-stiffened panel, ju y D x can be neglected when compared 
with D xy . In this case D 3 = D xy . The twisting stiffness D xy is calculated from Bredt 
theory (see, for example, ref. 24) and is given by 
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D 3 = D xy = 


2A 2 


bj_ + ^ + b 3 


(A24) 


{ a 66 x a 66 2 A 66 3 


in which A is the area enclosed by the hat and is given by 
A = + ^3) |p + ^4(0°) + t 2 (±45°)j 


(A25) 


Experimental work presented in reference 23 and analytical studies carried out to 
support reference 25 indicate that the value of Djjy given in equation (A24) is too large. 
In the studies presented in this report the twisting stiffness was obtained by applying a 
0.3 factor to the value given by equation (A24). 

In the case of the corrugated panel, PyD x cannot be ignored when compared with 
Dxy. The stiffness D3 is given by 

D = 2b 2( D 12 2 + 2D 66 2 ) + 2b 3( D 12 3 + 2D 66 3 ) 

3 p 
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TABLE I.- EQUATIONS FOR STRESS ANALYSIS 


Loading-equation 

identification 

Equation 

Equation 

number 

Longitudinal 

e = Nx 
x ET 

( 1 ) 

compression 

A 12 k 
' y >'‘ £x A 22k 

( 2 ) 


Nx k = e x ET k 

(3) 

Shear 

N xy 4 = N xy 

(4) 


Nxy 

Nxy 2 - N xy 3 - ^ i 2 b 2 A 66l b 3 A 66l 
b l ^62 b l A 66 3 

(5) 


N xyj = N xy ' N xy 2 

( 6 ) 


y _ Nxy k 

A 66k 

(7) 

Lamina strains 

€j = cos 2 9 e x + sin 2 9 e y + cos 9 sin 9 y xy 

(8) 

referred to 


lamina principal 

e 2 = sin 2 9 e x + cos 2 0e y - cos 9 sin 9 y 

(9) 

axes 

y 12 = -2 sin 9 cos 9 e x + 2 sin 9 cos 9 e y 



+ (cos 2 9 - sin 2 9 ) y 

xy 

( 10 ) 

Lamina stresses 

°\ = Qn e i + Ql 2 e 2 

(ID 

referred to 
lamina principal 

a 2 = Q 12 e l + Q 22 e 2 

( 12 ) 

axes 

t 12 = Q 66 V 12 

(13) 
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TABLE H.- EQUATIONS FOR STRENGTH FAILURE CRITERIA 


Failure 

criterion 

Equation 

(*) 

Equation 

number 

Maximum 

strain 

e l 

ih 1 

(14) 


e 2 

r* sl 

2 

(15) 


y _U sl 

y a 
12 

(16) 

Stress 

interaction 

a l 

■i * 1 

CT 1 

(17) 


/ X 2 /T \2 

fe si 

\ a 2/ \12/ 

(18) 


* 

Superscript a denotes allowable value. 



o> 


TABLE m.- EQUATIONS FOR OVERALL AND LOCAL BUCKLING 


Loading 

Equation 

Equation 

number 

References 

Overall buckling 

' Longitudinal 
compression 

N x crit = \ 

X ’ Cnt dL 2 j , »2EI 

2L 2 b 2 A 6 6 2 

(19) 

Equation (92’), ref. 11; 
equation (3), ref. 12 

itoiDo 

Shear Q = - 

; ° 3 

, | For 0 > 1, Nxy)Crit = 8.125 + 

. For 0 < 1, Nxy^rtt = + 0.5320+ O.9380 2 ) 

(20) Equations (2.2.2-21) 

and (2.2.2-22), 

: ref. 8; see, also 
(21 ^ j pp. 383-384, ref. 11; 

pp. 468-471, ref. 13 

(22) 

Combined 

loads 

Nv ! N xv \ 2 

I (23) j Equation (105.8), ref. 13; 

| i refs. 14 and 15 

N XjC rit V N xy,crit/ 

Local buckling; equations applied to each element 

Longitudinal 

compression 

N x ,crit = ^|l>llD22) 1/2 + »12 + 2D 66 


(24) 

Equation (2.2.2-16), ref. 8; 
equation (233), ref. 11 

Shear 

, ( D 11 D 22) 1/2 
d 12 + 2D 66 

For 0>1, = (|) 2 (D n D \ 2 f\ 

For 9 < 1, Nxy,crit ~ (j") 1^22 ( D 12 + 2 

8.125 + 

e 1 

D 6e)J 1 ^ 2 (l 1 - 7 + 0.5320 + O.9380 2 ) 

(25) 

(26) 
(27) 

Equations (2.2.2-21) 
and (2.2.2-22), 
ref. 8; see, also 
pp. 383-384, ref. 11; 
pp. 468-471, ref. 13 

Nx ( N xv \ 3 

Combined Loading critical when rj + 1— j = 1 (28) Equation (105.8), ref. 13; 

loads x > crit V N xy,crit/ ref . 15 


Combined 

loads 


Loading critical when 


(28) 


Equation (105.8), ref. 13; 
ref. 15 



TABLE IV.- PROPERTIES OF GRAPHITE -EPOXY MATERIAL 

USED IN SAMPLE CALCULATIONS 

[ten. denotes tension; comp, denotes compression; and 
superscript a indicates allowable] 


Symbol 

Value in 
SI units 

Value in U.S. 
Customary units 

Density and elastic properties 

P 

1520 kg/m3 

0.055 lbm/in3 

Ej 

145 GPa 

21 x 10 6 psi 

e 2 

16.5 GPa 

2.39 x 10 6 psi 

G 12 

4.48 GPa 

0.65 x 10 6 psi 

Pl2 

0.314 

0.314 

P 21 

0.037 

0.037 


Allowable stresses 

cr| (ten.) 

1240 MPa 

180 ksi 

cr| (comp.) 

-1240 MPa 

-180 ksi 

o| (ten.) 

55.2 MPa 

8 ksi 

(comp.) 

-207 MPa 

-30 ksi 

T a 

12 

82.7 MPa 

12 ksi 


Allowable strains 

(ten.) 

0.0087 

0.0087 

(comp.) 

-0.0087 

-0.0087 

e| (ten.) 

0.00475 

0.00475 

e| (comp.) 

-0.01764 

-0.01764 

y 3- 

r 12 

0.01846 

0.01846 



TABLE V.- PROPERTIES OF ALUMINUM USED 


IN SAMPLE CALCULATIONS 


Material properties 

Symbol 

Value in 
SI units 

Value in U.S. 
Customary units 

P 

2780 kg/m 3 

0.101 lbm/in 3 

E 

72.4 GPa 

10.5 x 10® psi 

G 

27.4 GPa 

3.98 x 10 6 psi 

P 

0.32 

0.32 

r a 

221 MPa 

32 ksi 



TABLE VI.- LIMITS ON DESIGN VARIABLES 


Dimension 

Hat 

-stiffened 

panel 

Corrugated 

panel 


cm 


in. 

cm 

in. 

t 3 (0°), maximum 

0.318 


0.125 

0.318 

0.125 

t 4 ( 0 °), maximum 

.635 


.250 

Not applicable 

minimum 

2.03 


.80 

2.03 

.80 

b 2 ? minimum 

1.02 


.40 

1.02 

.40 

b 3 , minimum 

2.03 


.80 

2.03 

.80 

b 4 , minimum 

1.02 


.40 

Not applicable 
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